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Abstract 

In the pyrochlore compounds, Tb2Ti207 and Tb2Sn2C>7, only the Tb 3+ ions are magnetic. They exhibit quite 
abnormal - and, in view of their chemical similarity, strikingly different - magnetic behaviour, as probed by neutron 
diffraction at ambient and applied pressure. Tb2Ti207 is a cooperative paramagnet ('spin liquid'), without long 
range order at ambient pressure; however, it does become ordered under pressure. By contrast, Tb2Sn207 enters 
an "ordered spin ice" state already at ambient pressure. We analyse a simple model which already clearly exhibits 
some of the qualitative features observed experimentally. Overall, comparing these two compounds emphasizes the 
power of small perturbations in selecting low-temperature states in geometrically frustrated systems. 
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1. Introduction 

Geometrical frustration (GF) [1] is now widely stud- 
ied in solid state physics, as it seems to play a key role 
in original phenomena recently observed in new mate- 
rials. Examples include the large anomalous Hall effect 
in ferromagnetic pyrochlores or spinels [2], the uncon- 
ventional superconductivity observed in water substi- 
tuted Naa;Co02 with triangular Co sheets [3], or the 
interaction between electric and magnetic properties 
of multiferroics materials [4] . 

What is geometrical frustration? Most simply, it oc- 
curs when the specific geometry of the lattice prevents 
magnetic interactions from being satisfied simulta- 
neously. In insulating systems such as the rare earth 
pyrochlores, the impossibility of a simple Neel ground 
state due to GF offers the possibility of finding a large 
variety of alternative, magnetic and non-magnetic, 
short- or long-ranged ordered states. In the most ex- 



* Corresponding Author Email : mirebca@dsm- 
mail.saclay.cea.fr 



treme case, paramagnetic behaviour persists down to 
the lowest temperatures, leading to an extended coop- 
erative paramagnetic, or spin liquid, regime, in which 
only short-range correlations result [5] . 

Ferromagnetic interactions on the pyrochlore lattice 
may also be frustrated, namely when the exchange 
is dominated by a strong anisotropy which forces the 
spins in a tetrahedron to point along their local, non- 
collinear easy axes [6]. This leads to the spin ice state, 
whose degeneracy can be mapped onto that of real ice 
[6], leading to approximately the same entropy in the 
ground state [7] . 

In real compounds, the eventual choice of the stable 
magnetic state depends on a subtle energy balance be- 
tween the frustrated first neighbour exchange energy 
term and perturbation terms of various origins (longer 
range interactions, anisotropies, quantum fluctuations, 
...). It is of course also determined by thermodynamic 
parameters, such as temperature, pressure or magnetic 
field. Counterintuitively, thermal fluctuations can even 
induce order when ordered states permit softer fluctua- 
tions than generic disordered ones. This effect is known 
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as order by disorder [8] and is commonly encountered 
in frustrated magnetism. It has been well studied by 
Monte-Carlo simulations, and also received some ex- 
perimental confirmation [9]. Pressure can change the 
nature of, and balance between different terms in the 
Hamiltonian, as they can depend on interatomic dis- 
tances in different ways. An applied field adds Zeeman 
energy, and can, for example, stabilize a subset of the 
original ground states, at times resulting in magneti- 
zation plateaus. 

In this paper, we study a well known pyrochlore 
Tb2Ti2 07, which we investigated by neutron diffrac- 
tion under extreme conditions of temperature (down 
to O.IK) and applied pressure (up to 8.7 GPa). We 
review one of its most fascinating properties, namely 
its ability to "crystallize" or order magnetically under 
pressure and we propose a new theoretical approach 
which accounts for some important peculiarities of this 
effect. We also compare Tb2Ti2C>7 to its sibling com- 
pound Tb2Sn2C>7, very recently studied, which behaves 
as an "ordered spin ice" . Both compounds have a fully 
chemically ordered structure, the pyrochlore structure 
of cubic Fd3m space group, where the Tb magnetic 
ions occupy a GF network of corner sharing tetrahe- 
dra. Although they differ only by the nature of the non- 
magnetic ion (Ti/Sn), they show very different mag- 
netic ground states. The comparison sheds some light 
on how to select the ground state through very small 
perturbations, one of the most prominent characteris- 
tics of geometrical frustration. 



2. Tb 2 Ti 2 07 : a spin liquid orders under 
applied pressure 

Tb 2 Ti2 07 is a famous example of a spin liquid, in- 
vestigated by numerous groups, where short range cor- 
related Tb spins fluctuate down to 70 mK at least [10], 
that is more than 300 times below the typical energy 
scale of the magnetic interactions (the Curie- Weiss con- 
stant 9cw of -19 K, where the minus sign corresponds 
to AF interactions). The persistence of these fluctua- 
tions was checked by muon relaxation [10], at the time 
scale of the muon probe of about 10~ 6 s. At shorter time 
scales, inelastic neutron scattering showed a quasi elas- 
tic signal, whose energy linewidth strongly decreases 
below about 1 K, indicating a stronger slowing down in 
this temperature range [11]. Coexisting with the spin 
liquid phase, spin glass like irreversibilities and anoma- 
lies of the specific heat were recently observed in the 
range 0.1 K-0.8K [12]. Using high pressure powder neu- 
tron diffraction [13], we observed two interesting phe- 
nomena induced by pressure [14] i) the onset of antifer- 
romagnetic long range order below a Neel temperature 
Tn of about 2 K: ii) the enhancement of the magnetic 



correlations in the spin liquid phase above Tn- Just be- 
low Tn , the ordered phase coexists with the spin liquid 
in a mixed solid-liquid phase, whose relative contribu- 
tions vary with pressure and temperature. The mag- 
netic Bragg peaks of the simple cubic lattice can be 
indexed from the crystal structure of Fd3m symme- 
try, taking a propagation vector k=(l,0,0). It means 
that in the cubic unit cell with four Tb tetrahedra, 
two tetrahedra are identical and two have reversed mo- 
ment directions. A longer wavelength modulation of 
this structure, involving a much larger unit cell, was 
also observed in the powder data. 

What is the pressure induced ground state? More 
fundamentally what is the role of pressure? To answer 
these questions, we performed single crystal neutron 
diffraction down to very low temperatures (0.14 K), 
combining hydrostatic pressure with anisotropic stress 
[15]. We showed that both components play a role in 
inducing the long range order, and that the ordered 
moment and Neel temperature can be tuned by the di- 
rection of the stress. A stress along a [110] axis, namely 
along the direction of the first neighbour distances be- 
tween Tb ions, is the most efficient in inducing mag- 
netic order (Fig. 1). 

FullProf refinements of the single crystal data al- 
lowed a determination of the magnetic structure with 
better precision, especially the local spin structure 
within a Tb tetrahedron. The structure corresponding 
to the best refinement (Rf=14% is given in table 1. 
The bond 1-4 along the axis of the stress, which should 
be reinforced, has AF collinear spins. This corresponds 
to a natural expectation for AF first neighbor ex- 
change. The orientation of the spin 2 (orthogonal to 
the 3 others) is more surprising since with 3 collinear 
spins, the exchange field on the fourth one should be 
also collinear. Since no collinear structure gave a good 
fit to the data, it suggests that the real spin struc- 
ture may be even more complex than the proposed 
one. In any case, both powder and single crystal data 
yield an important conclusion: we found that inside a 
tetrahedron, the magnetization is not compensated, 
namely the vectorial sum of the four spins is non zero, 
(although it is of course compensated within the cubic 



site x y z M 

1 0.5 0.5 0.5 [1 -1] 

2 0.25 0.25 0.5 [1 1] 

3 0.25 0.5 0.25 [-1 1] 

4 0.5 0.25 0.25 [-1 1] 

Table 1 

Orientation of the magnetic moments in one tetrahedron in the 
pressure induced state of Fig 1., deduced from the re incmcnt 
of the magnetic structure. The stress component is along [0 1 
1], The atomic coordinates x, y, z, are expressed in the cubic 
unit cell containing 4 tetrahedra. Two tetrahedra arc identical 
and two have reversed spin orientations. 
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Fig. 1. Tb2Ti207 : an antifcrromagnctic ordered state with 
k— (1,0,0) propagation vector is induced under pressure. Here 
an isotropic pressure Pi— 2.0 GPa is combined with uniaxial 
pressure P u —0.3 GPa along [011] axis. The variation of the 
peak intensity of the magnetic Bragg peak (120) shows the 
Neel temperature. The local spin structure in a tetrahedron 
has non compensated magnetization. 

cell, since magnetisations of the four tetrahedra cancel 
two by two). This means that in the pressure induced 
ground state, the local order does not correspond to 
any configuration which mimimizes the energy in the 
spin liquid phase. In other words, pressure does not 
select any energy state among those belonging to the 
ground state degeneracy of the spin liquid (the ground 
state expected if one considers Heisenberg spins cou- 
pled via first neighbour AF exchange interactions 
only). The anisotropic component (stress) relieves 
the frustration in a more drastic way, by creating 
uncompensated bonds, associated with a very small 
distortion of the pyrochlore lattice. In addition the 
isotropic component shortens all distances in the same 
way, increasing the frustrated exchange interaction. 
This effect could also contribute to the increase of Tn • 
The prominent role of stress in inducing magnetic 
order raises a subsequent question. Could it be sta- 
bilized spontaneously by internal stresses? To answer 
this question, we have now checked the magnetic or- 
der at ambient pressure by neutron diffraction in two 
Tb2Ti207 samples with different heat treatments, 
down to about 0.1 K. In an "as cast" powder sample, 
we observe at 0.07 K broad magnetic peaks close to the 
positions expected for the pressure induced magnetic 
order (Fig. 2). The Lorentzian lineshape corresponds 
to a finite correlation length of about 25 A (2-3 cubic 
cells). The peaks disappear around 0.3 K. We also 
studied a single crystal, which was annealed at 1150 
°C for 25 hours to relieve internal stresses. In the sec- 
ond case, the mesoscopic magnetic order is absent and 
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Fig. 2. Tb2Ti207 : a mesoscopic order is induced by spon- 
taneous strains at very low temperature. Magnetic neutron 
diffraction spectra at 0.07 K, showing broad peaks close to the 
positions of the (100) magnetic peak and secondary magnetic 
phase (S) of the pressure-induced state [14], The spectrum of 
the spin liquid regime at 1.2 K has been subtracted. The inci- 
dent neutron wavelength is 2.52 A. The broad magnetic peaks 
disappear at about 0.3 K. 

only the liquid-like correlations are observed, down to 
the minimum temperature of 0.15 K. Since both sam- 
ples are chemically ordered and stoichiometric within 
the accuracy of neutron diffraction, it means that the 
mesoscopic order is induced by internal stresses. The 
onset of this mesoscopic order may strongly influence 
the spin glass irreversibilities and anomalies of the spe- 
cific heat observed in the same temperature range[12], 
which seem to depend on the heat treatment. 



3. Tb 2 Sn 2 07: an ordered spin ice state 

In contrast to Tb2Ti207, Tb2Sn207 undergoes a 
transition to an ordered state already at ambient pres- 
sure. The magnetic structure, very recently determined 
by powder neutron diffraction experiments [16] was 
called an "ordered spin ice" . The local order within one 
tetrahedron is close to the "two in-two out" configu- 
ration of spin ice, taking into account a small devia- 
tion of 13°. of the magnetic moments with respect to 
the local <111> easy anisotropy axes. In the canoni- 
cal spin ice state, individual tetrahedra keep the mu- 
tual orientational disorder allowed by the "ice rules" , 
leading to short range order and ground state entropy 
[17]. Here the four tetrahedra of the unit cell are identi- 
cal, leading to an ordered structure with k=0 propaga- 
tion vector (Fig. 3). The resulting magnetic structure 
is non-collinear, but exhibits a ferromagnetic compo- 
nent, which represents about 37% of the Tb 3+ ordered 
moment. This explains the ferromagnetic character of 
the transition, previously observed by magnetization 
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Fig. 3. Tb2Sn207 : an ordered spin ice structure: the local 
spin structure in a tetrahedron is close to the "2 in-2 out" 
structure of a spin ice, but individual tetrahedra arc identical, 
leading to an ordered structure with k— propagation vector 
and ferromagnetic character. 

[18]. 

Together with the non-collinear magnetic structure, 
the original effects of the frustration persist in the or- 
dered phase of Tb2Sn207. The magnetic order is sta- 
bilized in two steps (1.3 K and 0.87 K) corresponding 
to anomalies of the specific heat, and not in a classi- 
cal second order transition. The correlation length in- 
creases throughout the transition region, and remains 
limited to 180 A even at very low temperature. The or- 
dered state coexists with slow collective fluctuations, 
in the time scale of 10 -4 -10 -5 s. They were probed by 
comparing the Tb 3+ moment value of 5.9(1) [j,b de- 
duced from neutron diffraction to the much lower value 
of 3.3(3) fiB deduced from the specific heat. 

The magnetic order in Tb2Sn207 may be compared 
to that found by Champion et al. [19], who consid- 
ered the competition between first neighbor exchange 
and uniaxial anisotropy in a pyrochlore ferromagnet. 
The model involves two parameters, the strength of 
the ferromagnetic interaction J and that of the uni- 
axial anisotropy D a along <111> axes. Ferromagnetic 
and spin ice states correspond to the cases D a /J = 
and D a / J = oo, respectively. For finite D a /J values, 
the magnetic order shows many similarities with the 
observed one. Namely: i) the ground state is ordered 
in a k=0 four sublattice structure, ii) the local order 
within one tetrahedron may also be deduced from the 
spin ice structure, iii) the magnetic transition is of first 
order, changing to second order with decreasing D a /J. 

However, the deviations from the local spin ice struc- 
ture are different in the model and in the real sys- 
tem. In the model, spins are uniformly canted towards 
the ferromagnetic direction. The ground state magne- 
tization relative to the local moment increases from 



0.578=1/73 (the average magnetization of a tetra- 
hedron in the spin ice case) to 1 (the ferromagnetic 
case) with decreasing D a /J. By contrast, in Tb2Sn207, 
the deviations of the magnetic moments from the lo- 
cal <111> axes actually reduce the magnetization (to 
about 0.37 in relative units). So the deviations of the 
magnetic moments from the local spin ice structure 
act in an opposite way to that predicted by the finite 
anisotropy ferromagnetic model. 

Finally, in Tb2Sn207, the neutron and magnetic 
data together with the comparison with theory, sug- 
gest that here the effective first neighbor interaction 
becomes ferromagnetic, although the physics of the 
system cannot be simply reduced to the energy scheme 
assumed in ref. [19]. 



4. A simple model for the influence of stress 

An analysis of the effect of pressure applied to an 
individual tetrahedron already manifests qualitatively 
two basic experimental results observed in Tb2Ti207: 
the much stronger influence of an anisotropic stress 
along a [110] axis (as opposed to hydrostatic pressure or 
to a stress along the [100] direction) , and the presence 
of an uncompensated magnetisation. 

In the isotropic problem, all six bonds of the tetrahe- 
dron are equivalent. Application of stress in the [110] 
direction lowers this symmetry, as shown in Fig. 4. 
In symmetry terms, the bonds form a six-dimensional 
representation of the tetrahedral group Td, which de- 
composes into three irreducible representations. A sin- 
glet, A amounts to a uniform change of all bonds to- 
gether. Furthermore, there is a doublet, E, the com- 
ponents of which correspond (a) to strengthening two 
opposite bonds, and weakening the four others (or vice 
versa) or (b) to weakening an opposite pair of those 
four bonds, and strengthening the other pair. Finally, 
each component of the triplet, T, implies a strength- 
ening/weakening (by an equal amount) of an opposite 
pair of bonds [20] . 

The crucial point is that the uniaxial [110] stress 
can couple to all three representations. In its pres- 
ence, there are three (instead of only one) symmetry- 
inequivalent bond strengths (see Fig. 4). In other 
words, the Hamiltonian including the uniaxial [110] 
stress has a lower symmetry than the isotropic one. 
The case of [1 0] pressure is intermediate: here 
only two representations are present, the triplet being 
absent, as illustrated in the left panel of Fig. 4. 

Whereas the initial degeneracy of the isotropic sys- 
tem is a signature of the different possible compromises 
between which bonds to frustrate and which to satisfy, 
some of these choices have become forbidden as it is 
not possible to trade off inequivalent bonds against one 
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Fig. 4. Tetrahedron under uniaxial stress, denoted by the ar- 
rows. Left (right) panel: stress applied in the [100] ([110]) 
direction. This splits the six bonds into the following sym- 
metry-incquivalcnt groups: the bond along the [110] direction 
(bottom), the one perpendicular to it (top) — which remain 
equivalent for [100] but not for [110] stress — and the four re- 
maining ones. 

another. 

For simplicity, let us consider a classical, isotropic 
Heisenberg antiferromagnet at T=0 under stress. We 
have considered the spin configurations which mini- 
mize the energy for different orientations of the stress. 
We find that these configurations have a compensated 
magnetization for a stress along [100]. By contrast, 
a non-compensated magnetic moment can arise for a 
stress along a [110] axis. A summary of the calculation 
is given below. 

The results for an ice-type model (i.e. a ferromag- 
net in the presence of anisotropy D a ) can be obtained 
along similar lines. It needs to be borne in mind, how- 
ever, that (a) the strict ice model (D a /J=oo) does not 
permit small deviations of the spins from their pre- 
ferred axes, and that (b) a ferromagnet will generically 
exhibit a non-compensated moment even in the ab- 
sence of stress. By contrast, an anisotropic (but strain- 
free) antiferromagnet has a momentless ground state, 
which, however, is the simple FeF3-type 'all-in' or 'all- 
out' ground state. 

This classical isotropic Hamiltonian has a continuous 
two-parameter family of degenerate ground states in 
the isotropic case [1]. This degeneracy is reduced by the 
strain. For example, in the presence of an E-distortion 
weakening the average strength of the top/bottom pair 
of bonds with respect to the other two pairs, a collinear 
state will be selected. In this state, each spin is parallel 
to its partner at the other end of the coloured bond, 
and antiparallel to the other pair of spins. The total 
spin of the tetrahedron thus remains compensated at 
zero. For an E distortion of the opposite sign, the top 
(and bottom) pair of spins will be antialigned; for the 
full pyrochlore lattice, this generates decoupled chain 
states. In contrast to this situation, the presence of a T 
distortion does not change the energy of the isotropic 
ground states relative to one another to first order. 
This happens because it couples to a difference in the 
expectation value of opposing bonds, a difference which 
vanishes in the unperturbed ground states. 

However, in higher order, a difference in relative 



bond strength can induce a difference between the ex- 
pectation values of the scalar product Si ■ Sj across the 
bottom and top bonds. Such a difference is equivalent 
to an uncompensated total moment of the tetrahedron 
(as a tetrahedron with zero moment necessarily has 
equal expectation values of Si ■ Sj on opposite bonds) . 
The ground states in the presence of stress are thus 
close to - but not a subset - of those of the isotropic 
system. 



5. Discussion 

In this section we briefly comment about the rele- 
vance of the above model to Tb2Ti2C>7, and then focus 
on the origin of the differences between the two com- 
pounds. 

At ambient pressure, the fact that Tb2Ti2C>7 does 
not order but remains 'liquid' down to 70 mK, is still a 
challenge to theory. Sophisticated calculations taking 
into account the crystal field energy [21] together with 
dipolar interactions, predict an Ising like behavior for 
the Tb ion moment in the the ground state (with the 
moment reduced with respect to the free ion value) and 
an effective AF first neighbor interaction [22]. These 
calculations predict at ambient pressure a transition 
to an AF order similar to that found in FeF3 (a k=0 
structure with an "all in-all out" local configuration) 
below 1-2K, which is however not what is observed 
experimentally. 

The simple model discussed above, already proposed 
on a more empirical basis in Ref. [16], reproduces the 
main characteristics of the pressure- induced state in 
Tb2Ti2C>7, namely the stronger effect of degeneracy 
lifting of the [110] over the [100] stress, and the appear- 
ance of an uncompensated magnetisation in the for- 
mer case. This is presumably the case because it incor- 
porates the most fundamental property of the stress, 
namely the explicit symmetry breaking it induces. This 
effect should occur in qualitatively the same way in a 
much larger class of models. 

These results are therefore rather robust but, by the 
same token, they are also only qualitative: the model 
in its current form yields little information on the de- 
tailed Hamiltonian of the system, nor the origin of the 
effective nearest-neighbour interaction J and of its sen- 
sitivity to pressure. In particular, we have not been 
able to reproduce the detailed finite-temperature spin 
structure. 

We now turn to Tb2Sn207, which behaves as an or- 
dered spin ice. Neutron data as compared with theory 
strongly suggest that the effective first neighbour in- 
teraction has now become ferromagnetic. What is the 
reason for this change? We can propose the follow- 
ing explanation. In the "true" spin ices (Ho2Ti2C>7 or 



5 



Dy2Ti2 07 with stronger uniaxial anisotropy), it was 
shown that the effective ferromagnetic interaction re- 
sults from the influence of the dipolar coupling which 
overcomes the weak AF superexchange [17] . Taking the 
same conventional notations, the effective first neigh- 
bour interaction J e // is expressed as J e //= Jnn+D„„, 
where J n „=J/3 and D nn =5D/3 are the superexchange 
and dipolar energy scales, respectively. In Tb2Ti207 
(Jnn=-0.88 K, D n „=0.8 K from ref. [22]), this effective 
interaction remains AF. In Tb2Sn207, Sn substitution 
enlarges the unit cell (from a= 10.149 to 10.426 A in 
Ti and Sn compounds respectively). This expansion 
Aa/a ~ 2.7%, equivalent to a negative chemical pres- 
sure of about 12-15 GPa, should strongly decrease the 
AF superexchange interaction J. Assuming roughly a 
decrease of J nn in the ratio of the Curie-Weiss con- 
stants (-19 K and -12 K in Ti and Sn compounds re- 
spectively) without big changes in the dipolar constant, 
we get J e // = 0.18 K > for Tb 2 Sn 2 7 . Therefore the 
expansion in the unit cell induced by Sn substitution 
might be enough to switch the compound to the spin 
ice region of the phase diagram [23] . 

To go further, microscopic models should take into 
account the exact nature of the anisotropy, which is not 
simply uniaxial in the Tb compounds [21]. This involves 
a reinvestigation of the crystal field levels, currently 
in progress. It could exhibit more subtle differences 
between the two compounds than the simple effect of 
a chemical pressure discussed above. 

In conclusion, the two compounds studied here 
clearly show the rich variety of behaviour exhibited by 
geometrically frustrated magnets. Comparing them 
allows one to understand better the key role played 
by small perturbations in selecting one peculiar state 
among the many potential magnetic states. 
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